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Abstract 

For a particle in orbit about a static spherically symmetric body, we study the change in self-force 



that results when the central body type (i.e., the choice of interior metric for the Schwarzschild 
exterior) is changed. While a straight self-force is difficult to compute because of the need for 
regularization, such a "self-force difference" may be computed directly from the mode functions 
of the relevant wave equations. This technique gives a simple probe of the (non)locality of the 
(3JT)! force, as well as offers the practical benefit of an easy determination of the self-force on a body 



orbiting an arbitrary (static spherically symmetric) central body, once the corresponding result for 
a black hole (or some other reference interior) is known. We derive a general expression for the 
self-force difference at the level of a mode-sum in the case of a (possibly non-minimally coupled) 
scalar charge and indicate the generalization to the electromagnetic and gravitational cases. We 
then consider specific choices of orbit and/or central body. Our main findings are: (1) For charges 
held static at a large distance from the central body, the self-force is independent of the central 
body type in the minimally coupled scalar case and the electromagnetic case (but dependent in 



X 

the nonminimally coupled scalar case); (2) For circular orbits about a thin-shell spacetime in the 
scalar case, the fractional change in self-force from a black hole spacetime is much larger for the 
radial (conservative) force than for the angular (dissipative) force; and (3) the radial self-force 
difference (between these spacetimes) agrees closely for a static charge and a circular orbit of the 
same radius. 
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I. INTRODUCTION 



The leading order force exerted on a small body by its own self-field is known as the 
self-force. Heuristically, this force can be broken up into two pieces. One is local: the body's 
motion (and/or the properties of the medium it moves through) may distort the local self- 
field, causing an asymmetry over the body that exerts a net force. The other is non-local: 
field "sent out" by the body may encounter the body again at a later time, resulting in 
a force that formally depends on the entire past history of the particle motion. This type 
of effect — called a "tail" effect — is unfamiliar because it does not occur for electromagnetic 
fields in empty four- dimensional flat spacetime, but it is fact generic: higher dimensional 
fields, massive fields, fields propagating on a curved spacetime, and fields propagating in the 
presence of non-trivial boundary conditions (such as caused by the presence of conductors) 
will generically give rise to non-local self-forces. 

The above heuristic split may be made precise through the use of the Hadamard decom- 
position (e.g., [lj]) of a Green's function, which locally defines "direct" and "tail" pieces, 
supported on and interior to (respectively) the future light cone of the source point. As 
demonstrated explicitly in the pioneering work of DeWitt and Brehme Q in the context of 
a charged particle moving in a fixed curved background spacetime, the direct part gives rise 
to a local self-force, while the tail piece — together with contributions from the region outside 
the normal neighborhood, where the decomposition is not defined — gives rise to an integral 
of the derivative of the Green's function over the past history of the charge. While we 
will consider one electromagnetic example in section IHIl in this paper we will be primarily 
concerned with the case of a scalar charge first analyzed by Quinn [3J], which shares many 
of the features of the electromagnetic case (and gravitational case) while avoiding much of 
the computational complexity. Quinn considers the motion of a point particle coupled to 
a scalar field on a fixed globally hyperbolic background spacetime. The spacetime is not 
required to be vacuum, but the use of the scalar wave equation (as opposed to some coupled 
matter and scalar field equation) entails the assumption that any matter present does not 
couple to the scalar field. The spacetimes of this paper will contain matter only away from 
the particle, in which case Quinn's force becomes 



- (a" - oV) + lim f (g^ + u^u v )W v G(z(r),z(T'))dT' 
3 e->o+ 
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where G(x,x') is the retarded Green's function for the scalar wave equation, satisfying 

[g^y^u ~ £R] G{x, x') = -4vr5 4 (x, x') (2) 

for some coupling constant £. While Quinn did not include this coupling to curvature, his 
local analysis extends trivially to £ 7^ when the particle moves through a vacuum region, 
as assumed here. We include the curvature coupling in order to make equation (j2J) more 
analogous to Maxwell's equation and the linearized Einstein equation, both of which involve 
couplings to curvature when expressed in wave-equation form (i.e., in the Lorenz gauge; 
c.f. (J3"T]) and e.g. [lj). In equation ([T]) we take an explicitly perturbative viewpoint, where 
the particle's lowest-order motion z^(t) (four-velocity w M , four-acceleration a M , and four-jerk 
d M ) is caused by an external scalar field (or other force), and the small self- force correction 
is viewed as a vector field defined on z M (r), to be used perturbatively or incorporated 
self-consistently in the spirit of 4, l5f| . Note that since mass is not conserved for scalar 

nn 

charges, there is also a self-field effect on the mass [3|, |6j. Since we view the scalar field 
primarily as a toy model for understanding the more physically interesting electromagnetic 
and gravitational cases (where mass is conserved), we do not compute this effect. Similarly, 
we assume that the scalar charge q is constant, although this is not required by the scalar 
theory Q, B, B] - 

The first term in equation ([T|) is the local self-force coming from the direct part of the 
Green's function. It takes an identical form to the classic Abraham-Lorentz-Dirac equation 
(c.f. ( |30|) and e.g. Q|) for electromagnetic fields, depending on the derivative of the acceler- 



ation. The second term is the nonlocal self-force coming from the tail portion of the Green's 

n 

function!!! Although this "tail integral" depends upon the entire past history of the particle, 
in physically reasonable situations one expects the contribution from the distant past to be 
negligible. That is, there must be some point in the past at which the tail integral can be 
cut off, with numerically negligible error. The question of just how far in the past one must 
integrate is the question of the locality of the tail integral. 

Several authors have already investigated this question js-12], primarily in the context of 
attempting to compute the tail integral via a series expansion of the Green's function in the 
normal neighborhood of the charge. A generic finding of these authors is that the region of 



1 The limiting procedure ensures that the direct portion of the Green's function — which is nonintegrably 
singular — is not included. See (e.g.) Poisson |l[ for details. 
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non-negligible support of the tail integral extends well beyond the normal neighborhood: the 
tail integral is too non-local to be handled by series expansions alone. This non-locality can 
also be seen in the case of a slow-moving particle in a spherically symmetric weakly-curved 
spacetime, in an approach originally due to deWitt and deWitt [13] in the electromagnetic 
case and later improved and extended to the scalar and gravitational cases by Pfenning and 
Poisson [14]. In these remarkable calculations, one considers a body that is freely falling 
to lowest order, so that the direct (local) portion of the self-force vanishes and only the 
tail (non-local) portion remains. However, since by assumption this body moves slowly 
in a weakly curved spacetime, one should be able to regard the effects of curvature as a 
"force" causing a gravitational "acceleration" g defined relative to inertial motion in the 
flat background. From this point of view one would expect the self-force to be given by the 
usual Abraham-Lorentz-Dirac type direct force, proportional to the time derivative of the 
"acceleration" g. Indeed, this result is recovered: when the tail integral is computed, it gives 
rise to exactly the Newtonian expectations. 

However, this result comes with a twist: the tail integral that gives rise to the local 
Newtonian answer is in fact highly non-local! Rather than it becoming dominated by a local 
contribution in the weak- field limit, as one might expect, it turns out that one must integrate 
at least a light reflection time (off the central body) in the past to recover the simple local 
self- force. In fact, if one views the central body as a point mass, then in the scalar case as 
clearly explained by Anderson and Wiseman [l2] — the tail integrand actually vanishes in 
the immediate past of the particle, giving its first contribution at a light reflection time off 
the point mass in the past. This provides a dramatic illustration of the way in which a local 
expression for the force can arise from a highly non-local tail integral. We will refer to such 
a self-force as weakly local: it depends only on local properties, but the tail integral may be 
non-local. 

A direct evaluation of the tail integral is not the only way to calculate a self-force. Since 
the tail integral appears as part of an analytic expression for the field near a point particle, 
it is often (and even usually) more convenient to simply compute the field of a point particle 
by some other method and then extract the tail piece by a regularization procedure. This 
type of approach gives no direct information on the locality of the integral, but by comparing 
the results of different calculations one can glean some relevant information. For example, 
consider two spacetimes that agree in a region of the charge but disagree elsewhere. If there 





star SF 


black hole SF 


scalar case 


F r = 2£q 2 M/r 3 





EM case 


F r = e 2 M/r 3 


F r = e 2 M/r 3 



FIG. 1. A table of known self- force results for static charges in the far field of Newtonian stars and 
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black holes. The Newtonian star results are taken from Pfenning and Poisson [14], while the black 



hole results (which are exact if r is the Schwarzschild coordinate) are due to Wiseman 



scalar case and Smith and Will 



161 ] in the 



151 ] the electromagnetic case. Here M is the mass of the central 



body, and q/e is the scalar /electric charge of the particle. 

is a disagreement in the self-force on a particle moving through the region of agreement 



then the tail must be at least as non-local as the light reflection time across the region 
And if the results agree for a large class of locally-agreeing spacetimes, one would say that 
the self-force is weakly local in that case. 

An example of this type of comparison was discussed by Pfenning and Poisson [ijj], 
who compared their results for Newtonian central bodies with known results for black holes 



m 



16]. This comparison is conveniently summarized in tabled which shows results for static 

charges in the far-field of Newtonian stars and black holes. We see that in the electromagnetic 

and minimally-coupled (£ = 0) scalar cases, the far-field self-force is independent of the 

central object type (Newtonian star or black hole), whereas for the nonminimally-coupled 

scalar case there is a difference. This difference confirms the non-locality of the tail integral; 

but it is interesting that weak locality holds only in the electromagnetic and minimally 

coupled scalar cases. What is so different about the nonminimally coupled scalar case? 

In order to better understand as well as significantly extend these types of results, in this 

paper we make a direct attack on the dependence of the self-force on the central object. In 

effect, we take the method of comparison to an extreme: rather than computing a self-force, 

we will directly compute the difference in self-force that results from changing the interior 

metric while keeping the Schwarzschild exterior the same. Such a "self-force difference" 

is easy to compute because the direct portions of the Green's functions cancel out upon 

2 More precisely, consider the causal diamond defined by the current particle position z(t) and a prior 

particle position z(t'). If this diamond lies within the region of agreement, then the Green's function 

G(z(t), z(t')) must agree for both spacetimes. Thus, any difference in the self-force must come from 

times t' far enough in the past that the causal diamond intersects the central body. 
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subtraction, automatically picking out the relevant tail portion. More precisely, for points 
x and x' within the region where the two spacetimes agree, define the difference AG(x, x') 
between the retarded Green's functions. But initial value results imply that the Green's 
functions agree for points whose causal diamond lies within the region of agreement, so that 
in particular the difference AG must vanish in a neighborhood of coincidence, x = x' . This 
means that when computing the difference AF M in force on a particle moving through the 
region of agreement, we may actually drop the limit in equation (JTJ, giving simply 

AF M = q 2 J {jf + uV) V u AG(z(r), z(r'))dr'. (3) 

Note that we have extended the worldline integral in the future direction, which is possible 
because of the vanishing of the support of the retarded Green's function for x' in the future 
of x. As we shall see, equation ([3]) can be evaluated directly from the mode functions of 
the relevant spacetimes, avoiding the delicate regularization needed to compute a straight 
self-force, equation (jTJ). 

We consider static spherically symmetric bodies (described by a metric of the form (j3D) 
and obtain the following results. First, we "fill in the gaps" between Newtonian stars and 
black holes in table [J providing analogous results for arbitrary central bodies. Specifically, 
we show that the far-field self-force on a static charge is indeed independent of the central 
body (i.e., weakly local) in the electromagnetic and minimally-coupled scalar cases, and that 
the force indeed depends in detail on the body type in the nonminimally coupled scalar case. 
In this latter case we provide a specific example of the detailed dependence by considering a 
central body that consists of a thin-shell of radius tq. As one changes the radius r$ between 
some very large radius (Newtonian object) and tq = 2M (black hole), one can see the 
far- field force interpolatqj between the two cases (figure [2]). Our analysis also sheds some 
light on the reason the nonminimally coupled scalar case behaves differently from the other 
cases: weak locality in this limit is identified with the property that the static spherically 
symmetric solution to a wave equation is the same for all metrics. This property holds 
for the minimally-coupled scalar and electromagnetic equations (the solution is simply the 
constant solution), but does not hold for the non- minimally coupled scalar wave equation. 

We next move to the simplest dynamic case of circular orbits. Here there is a dissipative 
self-force — the result of energy being radiated away from the system — in addition to the 

However, the transition is not smooth everywhere, because instabilities can occur in some regions of 

parameter space — see appendix. 



conservative effects seen in the static case. We numerically compute the self-force difference 
between a thin-shell spacetime and Schwarzschild (figure [3]), finding that the difference in 
interior is far more important for the conservative part of the force than for the dissipative 
part (figure H]). This is in accord with the intuition that the dissipative self-force is more 
like a local radiation reaction force, whereas the conservative self-force should be sensitive 
to boundary condition J This constitutes numerical evidence for weak .oca.it y of the di 8 - 
sipative force even in the strong-field regime, and holds out the possibility of a reasonable 
degree of true tail locality in this regime|f| We also notice (figure [5]) that the conservative 
self-force difference closely matches the static case self-force difference (whereas the straight 
self-forces do not match). 

Although this paper gives little quantitative handle on the precise degree of (non) locality 
of the tail, we take a moment to point out a potential practical application, should such 
knowledge ever be gained. Much current effort is directed towards t he p roduction of wave- 
form templates for extreme mass-ratio inspiral to be used in LISA |l7| data analysis, and 
here the gravitational self-force must take a central role (see 18j for a recent review). The use 
of frequency- domain techniques — convenient for circular orbits — becomes much less feasible 
for very complicated orbits, where a great many frequencies would be needed to reconstruct 
the field. However, suppose one knew that the tail integral was precisely T local, in that 
one could cut off the integration at a time T in the past with negligible error. This means 
that in an indirect computation of the force, one can replace the particle orbit with anything 
one wants in the region before time T in the past, while not affecting the value of the self- 
force. This freedom could potentially be exploited to greatly improve the convergence of the 
Fourier series used to reconstruct the field and force. For example, a zoom-whirl orbit of Kerr 
could be replaced by a (non-geodesic) circular orbit of the whirl radius and frequency. More 
generally, one could imagine having an algorithm to modify the past properties of orbits in 
order to maximize the estimated speed of convergence of the inverse fourier transform. 

Finally, we note that the techniques used in this paper could be used to help generate 
waveform templates for the case where the central body is not a black hole. For spinning 
objects, the exterior spacetime metric will depend on the central body, and any departure 

4 A simple example of a conservative self-force having nothing to do with curved spacetime is the force on 

a charged particle in the presence of a conductor. 

5 The idea that the tail may become more local in the strong-field regime can be supported by the intuition 

that field is bent more strongly (and hence quickly) back to the particle in the presence of large curvatures. 
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from a Kerr black hole will show up at lowest order in the motion of the inspiraling body. 
But for non-spinning bodies (that are also roughly spherical), the exterior metric will remain 
Schwarzschild and differences in motion will be seen first in the self-force. However, by the 
methods used in this paper, this self-force is easily computed for a given body type once 
the corresponding result for Schwarzschild is known. Therefore, it appears that it should 
be a relatively computationally inexpensive matter to broaden search templates to include 
central objects that are not black holes, even in the non-spinning case. 

In section [III we derive a general frequency-domain expression for scalar self-force dif- 
ferences involving static, spherically symmetric central bodies. In section III II we evaluate 
this expression (and its electromagnetic analog) for static charges in the far-field limit and 
extend the results of table [I] to arbitrary central bodies. Finally in section IIVI we consider a 
thin-shell central body and compute the self-force difference for static and circular orbits of 



arbitrary radius. Our conventions are those of Wald 
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II. GENERAL EXPRESSION FOR SPHERICALLY SYMMETRIC CENTRAL 
BODIES 

To construct a metric representing a general static spherically symmetric body we will 
match a general interior metric form (see e.g. [20]) with the Schwarzschild exterior. Taking 
the matching surface to be at a radius ro > 2M, our spacetime will be 

| -e^^dt 2 + e 2A ^dr 2 + fW < f < r 
ds 2 = , (4) 

I -f{r)dt 2 + /(r)-Mr 2 + r 2 dtt 2 r > r 
where f(r) = 1 — 2M/r and dQ 2 = d6 2 + sin 2 9d(f> 2 , and \l/ and A are analytiqj functions 
such that e 2 **- 1 " -* = /(r ). The last requirement arises from the demand that the induced 
metric agree on the matching surface (so that the geometry is smooth enough to interpret 



as a single solution 



21 



22] ), together with our choice to use the "same" coordinate t on 
both sides of the matching surface. However, this places no physical restrictions on the 
interior metric, since a rescaling of the time coordinate can always be used to rescale \1/ so 
as to satisfy e 2 *^ -* = /(ro). We do not demand that the extrinsic curvature matches, which 



corresponds to allowing a surface layer of matter to be present at r = r 



21 



22]. 



6 The assumption of analyticity enables us to establish existence and uniqueness properties for the mode 
functions associated with the interior metric. 
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While the metric components normal to the surface play no role in the matching analysis, 
it is often convenient to work with the same radial coordinate on both sides. In this paper 
we will make the simple choice that 

f = ro + ( r _ ro )e- A ^)/(r )- 1 / 2 , (5) 

given which the metric (jlj) may be written in the single coordinate system (t,r,8,<f)) as 

- e mf)dt 2 + e^^-^^fir^dr 2 + fW r [l - e A ^f(r ) 1/2 ] < r < r 
-f(r)dt 2 + /(r)" 1 ^ 2 + r 2 dtt 2 r > r 



ds 2 



(6) 

In these coordinates, all metric components agree at the matching surface r = r$. 

We will solve for the Green's function (J2]) on the background (j4]) by making use of sepa- 
ration of variables. The time translation and spherical symmetries of the metric imply that 
we should decompose as followsjll 

/oo oo £ 
duJ E E Y ^ WtJP, 4f)GUr, r'K^e-*', (7) 
'°° £=0 m=-£ 

where Y/ m are the spherical harmonics [23]. Note that we write Gi u since the rotational 
symmetry also implies that the mode functions are independent of m (as seen explicitly 
in equations (|SJ) and fl9]) below). After likewise decomposing the delta-function source into 
modes, equation (jSJ) implies 

f < r : G'i{f) + (-A' + I + GL(r) - e 2A - e~™u 2 + £r) G^{f) = 

(8) 

- > : Gl(r) + 2 -^pGl(r) - (^±11 _ ^ G&j(r) = _ilL 5(r _ r ^ (9) 

where we have defined r q = r' to allow the prime to indicate ordinary differentiation. (The 
notation r q reminds that the particle will be at the radius r = r q , which we assume to be at 
r > ro.) The solution for is constructed by matching together homogeneous solutions 
from the regions < f < ro, ro < r < r q , and r > r q . To determine the matching conditions 



7 There is no guarantee that this decomposition will always exist for the class of spacetimes fU we consider. 
For example, the presence of instabilities could render the mode functions G&j ill-defined, as happens for 
the case discussed in the appendix. We simply assume that this decomposition exists, in the hopes that 
this eliminates only unphysical spacetimes. 
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it is easiest to work with the same coordinate r on both sides of ro, i.e., to use the metric 
form 06]). In this case the wave equation can be written in the form 

Air 

G'i{r) + M(r)GL(r) + [N(r) + V(r - r )] G iu) (r) = —5(r - r g ), (10) 

where N(r) is a continuous function, M(r) is a continuous function except for a possible 
jump discontinuity at r , and x is the coefficient of the delta-function contribution to the 
Ricci scalar (coming from a possible boundary layer of matter at ro), 

rr +e 

X = lim / R(r)dr. (11) 

£ - >0 Jr -e 

It is now easy to see that the matching conditions are simply 

lim GUr) ~ lim G iu} (r) = (12) 

1 — Hq r— >r 

lim GL(r) - hm G^(r) = /"^J (13) 
lim G^(r) - lim G iu} (r) = (14) 

47T 



hmGL(r)- lim G'^r 



r2/ • (15) 

' r=r q 

Note that all derivatives are taken with respect to r (rather than f) in the above equations. 
We further impose the boundary conditions of regularity at f = and outgoing radiation 
{G^ — > r^ 1 e luir ) at r — > oo, which should correspond to choosing the retarded Green's 
function, as demanded by the self-force prescription ([T|). 

We will begin by performing the matching at r = r q , which is well-studied since equation 
([9]) is simply the Schwarzschild radial equation. In treating this equation it is convenient to 
choose homogeneous solutions Rf u and Rf° defined by the end behavioc 

RfJr) -> for r* ->• -oo (16) 

r 

RZ( r ) -> for r * -> +°°' ( 17 ) 

where r* is the "tortoise coordinate" defined by dr^/dr = / _1 and r* — > r as r — > oo. From 
equation ([7]) we see that R H corresponds to purely ingoing radiation at the horizon^ while 
8 The forms given in equations (116[) and (|17l) are not valid for oj = 0. In this case, one simply demands 
finite limits, so that Rf is regular at the horizon, and falls off at infinity. In fact, Rf Q and R^ are 



simply Legendre functions — see section HVl 
9 Although the Schwarzschild metric does not extend to r = 2M in our spacetime, it is still convenient to 

choose homogeneous basis solutions to © with properties at r — 2M in mind. 
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R°° corresponds to purely outgoing radiation at infinity. Matching these solutions at r = r q 
gives a useful particular solution to (jUJ), 

Att 

GST = -N^RUr<)RZ{r>l = r 2 fW{R H tR »y ^ 

where r < /r > indicates the lesser /greater of r and r q and W(ui,u 2 ) = U\u' 2 — u 2 u[ is the 
Wronskian of u\ and u 2 . Note that the differential equation ([9]) guarantees that the combi- 
nation r 2 fW(ui, u 2 ) is constant for any two homogeneous solutions u\ and u 2 , so that is 
a constant normalization factory] This particular solution to (Q has no incoming radiation 
either from infinity or the horizon, so that it should correspond to the retarded Green's 
function in pure Schwarzschild spacetime. This property makes it convenient for use in a 
computation of a self-force difference between our spacetime (HJ and the pure Schwarzschild 
spacetime. 

We now proceed to the matching at the surface of the body, r = r$. The general 
(inhomogeneous) solution to equation fl9]) compatible with our boundary conditions may 
be written as the particular solution ffT8j) plus a piece purely outgoing at infinity, 

GZ\r) = G$?"(r) + C £oj RZ(r), (19) 

where is a constant. On the other hand, the general solution to equation (jSJ) compatible 
with our boundary conditions is 

Gt(r) = DMr), (20) 

where D euj is a constant, and is the uniqui"1 (up to normalization) solution to equation 
(IE]) regular at the origin. The matching conditions ( |T2l) and ( |T3l) determine the constants 
Ce w and Out interest is in C^, which we determine to be 



n at dm/, \ Wjhu,, RZ) ~ W l XhJtfL 
Geoj = JytuKtofa)^^ 5^ t ;. u , r poo 



(21) 

r=ro 



In writing this expression we have assumed that the denominator in the fraction is non- 
zero. If the denominator does vanish, then our matching conditions (which follow from the 
equation for G&j) become inconsistent. Since the retarded Green's function must exist in 

that the Wronskian of R H and R°° is everywhere non-zero, so that this constant 



10 



It is also known (e.g. 



24 



is well-defined. 

11 The Frobenius method (see, e.g., 25j ) applied to equation (J8j) at the regular singular point r = implies 
that a unique (up to normalization) regular solution exists for ^> and A analytic at zero. 
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our globally hyperbolic spacetime 01]), this indicates that the mode-sum decomposition (J7J) 
is not defined in that spacetime. We show in the appendix that the denominator can vanish 
only for u = 0, in which case it corresponds to a linear instability of the spacetime. In the 
Wronskians present here and below, all derivatives are to be taken with respect to r (rather 
than f). 

Now consider the Green's function difference AG (defined only for r q > ro) between 
the general metric fll]) and pure Schwarzschild spacetime. Equation f|T9|) shows that this 
difference is controlled by the constant Ce u , so that we can write 

= N eu KUr )R%(r)RZ(r q ), (22) 

where 



wii^RD-u^xhvRl 



w{h UJ ,RZ)-U- l xh.RZ 



(23) 



r=ro 



is a constant (depending only on r ). Equations ( 122]) and f )23|) are the main results of this 
section, giving the "Green's function difference" AG mode by mode. Integrating the Green's 
function difference over a given orbit via equation ((3]) gives the difference in self-force between 
a body orbiting a static spherically symmetric body (equation fll])) and the same body on 
the same orbit about a Schwarzschild black hole of the same mass. As is made manifest in 
these results, this quantity can be computed directly from the mode functions I^.Rf^.RZ 
of the spacetimes, without the need for regularization. As previously discussed below (I2T|) . 
our expression for Kg^ is well-defined in all spacetimes (of the form fll])) for which the mode 
decomposition ([7]) is defined. 



A. Expression for two non-black-hole central bodies 

Suppose one is interested is the self-force difference for two different central bodies A and 
B (rather than for a body and a black hole). The most straightforward way to compute this 
difference is simply to use equations (122]) and fl23]) to obtain the difference from Schwarzschild 
for each central body, and then subtract. However, a more illuminating expression for the 
difference between A and B can be obtained as follows. Define Ae w to be the solution 
for central body A that is regular at the origin but now extended throughout the whole 
spacetime (j3J) by matching at r = r Q . Thus A £ul agrees with interior to the body and 
becomes a particular combination of Rf w and Rf° exterior to the body. By exactly the same 

12 



computation that leads to equation (1TB]) . the solution Gf w regular at the origin and outgoing 
at infinity is written (for r > tq) as 

G '- = - 4 V/u-'(^,i©' (24) 

where again the denominator is known to be constant by the properties of ([9]), and vanishes 
only in the presence of instabilities (see appendix and note at the end of this subsection). 
Defining the analogous quantities and Gf w for central body B, we arrive at the analogous 
expression. The difference A — B is then given by 

Since r 2 f'W (A^, R^) is known to be constant, we may choose any radius r to evaluate 
it (and likewise for B). In particular, we may choose r = r<, which enables us to use 
the property of the Wronskian that uiW{u2,u^) — U2W{ui,u^) = u^W{u2, u\) to simplify 
equation (1251) to 

\riAB atAB d oo /\ D oo / \ atAB a „ W(Be w ,Ae UJ ) . . 

^G lu) = -N iuJ R euj {r)R fMJ {r q ), N lul = 4tt— — — -— , (26) 

r J w K A iu>, R e J W [Ufa, R t J 

where N^f is a constant. This form shows that the difference in self-force induced by central 
bodies A and B is controlled by the Wronskian of their origin-regular mode functions A and 
B. In particular, when central bodies A and B are identical, then the mode functions A^ 
and Bfa will be linearly dependent, and the self-force difference vanishes by the Wronksian's 
vanishing on linearly-dependent functions. Equation (1261) also has the aesthetic advantage 
that it makes no direct reference to the "horizon" Schwarzschild solution Rf u , which should 
not play a privileged role in a self-force difference between two non-black-hole spacetimes. 
However, this formula offers no practical benefit over simply computing the difference from 
pure Schwarzschild by ( 1231) for each spacetime A and B, since the relationship of A^ (or 
Btu) to the corresponding solutions 1^ to the interior wave equation restores the complicated 
dependence on Rf u (or on some other choice of a second linearly independent solution to 
the Schwarzschild wave equation). Specifically, one has 

Ae.uj = E iuj Rfa + FfaRfa, (27) 
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where 



E 



W(I, 



oo 



W(Rg,R 



CO > 



r=ro 



w{i^Ri)-r^xh.R 



(28) 
(29) 



r=ro 



W(Rg,R?J 

Notice the simple relationship to the constant (equation fl23|) ) controlling the self-force 
difference between A and Schwarzschild; we have = —F^/E^. In particular equation 
is ill-defined when E iuJ = 0, which by (j2"7|) indicates linear dependence of A iuJ and R^. 



B. Electromagnetic and Gravitational Cases 

We now discuss the generalization of the above treatment of self-force differences to the 
electromagnetic and gravitational cases. The electromagnetic case proceeds in complete 
parallel with the scalar case. The electromagnetic analogs of ([T|) and (T5]) are 



= e 2 

and 



aV)+lim/ u u V lu G^\z(T),z(T'))u^(r')dT'\ (30) 

e ^ 0+ J-oo J 



[<*y V M V, - ir p ] G p a ,(x, x') = -A7cS a a ,S^(x, x'), (31) 

where the retarded Green's function should be chosen for G , J 12 l These equations give the 
self-force on a particle of charge e moving through a vacuum region of spacetime. As in the 
scalar case, consideration of a self-force difference allows one to drop the limit in equation 
( 130|) . producing the analog of equation ()3]). One can then follow the procedure used in 
section [Till to derive the analogs of (1711291). However, the separation of the equations is 
considerably more complicated than the scalar case, since (e.g.) vector spherical harmonics 
must be used to avoid mode-mode coupling. Since we perform no explicit calculations in the 
time-dependent case, we do not display these equations. Instead, in the following section we 
will develop the static case separately, where scalar harmonics can be used and the equations 
are considerably simpler. 

As in the scalar case, the electromagnetic self-force formalism is based on the wave equa- 
tion on a fixed curved spacetime, so that if the spacetime is non-vacuum, the results only 



12 Note that a vector potential constructed from G a a , using a conserved current will satisfy the Lorenz gauge 
condition, V^A' 1 = 0. 
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apply if the matter is not coupled to the electromagnetic field. While in the scalar case it 
seems fair to simply declare any matter to be uncoupled at a fundamental level, in the elec- 
tromagnetic case this would eliminate most known types of matter. Instead, the restriction 
on the type of matter in the electromagnetic case is best stated as the requirement of a 
nonpolarizable medium with no net charge. Thus, our electromagnetic results are valid only 
for nonpolarizable central bodies (with no net charge). The inclusion of polarizability is not 
difficult at a theoretical level — one would simply replace equation ( l3Tj) with a "macroscopic" 
Maxwell equation or other choice of coupled Maxwell and matter equations — but obtaining 
solutions for such equations may be more difficultcj We note, however, that in this case 
the self-force formalism would automatically include backreaction effects due to polarization 
induced by the charge; there is no clean distinction between these and other self-force effects. 
The theory of gravitational self- force [l[ B, 26-28] is also highly analogous to the scalar 



and electromagnetic cases, with the Linearized Einstein equation playing the role of the 
scalar wave equation / Maxwell's equation, and the perturbed geodesic equation "force" 
playing the role of the scalar force / Lorentz force. The main difference in the gravitational 
case is that the background spacetime is assumed to be vacuum. This is done because the 
assumption of having matter that does not couple to the field is no longer permissible, since 
gravity couples to everything. As explained and developed in [fjj], gravitational self- forces 
in the presence of matter require one to consider (linearized) coupled matter and Einstein 
equations. Thus, while the theory of gravitational self-force differences for static spherically 
symmetric bodies can be developed in analogy with the scalar and electromagnetic cases, it 
would require the use of Green's functions for (linearized) coupled matter and gravitational 
equations. 

III. ARBITRARY INTERIOR: STATIC CASE 

The case of a charge held static at a fixed radius r = r q outside the body provides 
enough simplification to obtain some general results (valid for a general interior metric). 



The worldline of the charge is given by u a = (\//(r 9 ), 0, 0, 0), whence equation ([3]) gives the 
13 An exception is the case of a conducting central body, where the matter model consists in the demand of 
zero electric field inside the body, which effectively contributes only a boundary condition to Maxwell's 



equation. Indeed, this case was considered by Shankar and Whiting 

Schwarzschild star), using a technique similar to ours. 
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29( (who also consider an insulating 



self- force difference to be 



AF r = q 2 ( v / 7^AG(r,r g ,a; = 0) 



(32) 



with all other components vanishing. However, Wiseman 16] has shown that the self-force 



on a static scalar charge outside a Schwarzschild black hole vanishes. Therefore, in the 
static case the self-force difference from Schwarzschild AF r is in fact equal to the self-force 
F r (and we will drop the A below). Notice that only the u = mode of the Green's 
function contributes to the static self-force. Accordingly, u> = will be implicit for the 
remainder of this section, and we will drop the label u on the mode functions. For oj = the 
Schwarzschild mode functions are known analytically 16|, |30] ; they are (up to normalization) 



simply the Legendre functions Pg(r) and Qe{r) [3l| of the Schwarzschild harmonic coordinate 
r = r /M - 1. Using the fact that W (P e (x) , Q e (x)) = 1/(1 - x 2 ) } equations (J22) and ([23]) 
become 

47T 

AG, = -—K e (r )Q e (f q )Q e (f) (33) 

and 



W(I t ,P l )-U- 1 xhPt 



(34) 



r=r 



w(i t ,Qi)-u-i x itQi 

where in equation ( 1341) and below it is understood that the Legendre functions are evaluated 
at f(r) (but that their derivatives in the Wronskian are taken with respect to r). Using 
equations © and ( 132|) . the self-force (difference) is seen to be 

,2 



AF r 



Q 



oo 

^~)Y^(2£ + l)K e (r Q )Q e (f g )Q' e (f q ), 



M 2 



(35) 



where the prime indicates differentiation with respect to r, and we have used the addition 



theorem 



23J to perform the sum over m. Equations (135]) and (|34|) give the self-force on a 



scalar charge held static outside an arbitrary central body of the form (J4j) . 

Now consider the "far-field" limit in which the charge is taken far away from the body, 
r q — > oo. For large x, the Legendre function Qe(x) behaves as so that the £ = 

mode dominates the far-field self-force. Specifically, we have 



q 2 M 



f(r q 



wii^p^-u^xhPv 



W(I ,Q )-£f-i X IoQo 



+ o 



r=r 



l M 3 



(36) 



We now prove that the first term in equation (1361) vanishes when £ = 0. Since Pq is simply 
a constant (non-zero) function, its Wronskian with I vanishes if and only if J is a constant 
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function. Now, Jo is the unique (up to normalization) solution to equation (JSJ) with i = u> = 
that is regular at the origin. When £ = 0, all terms proportonal to Gi u in equation (jSJ) vanish 
and the equation admits a constant solution, showing that Iq is constant and the Wronskian 
W(Iq, Pq) vanishes. Since the second term in the numerator of (|36|) also vanishes when 
£ = 0, the leading order far-field self- force vanishes when £ = 0. Thus we have the following 
result: 

• For minimally coupled fields, the leading order far- field self-force on a scalar charge held 
static outside an arbitrary static, spherically symmetric central body is 0(q 2 M 3 /r^). 
By contrast the force for non-minimally coupled fields is 0(q 2 M/r^). 

The result in the minimally-coupled case may be viewed as stating that the self-force is 
independent of the central body (always taking the value of zero) at 0(q 2 M/r^), confirming 
the suggestion of figure [H Similarly, we have shown that the corresponding force in the 
nonminimally-coupled case is dependent in detail on the central body composition. We note 
that some similar results have been obtained in [32j for expectation values of products of 
quantized scalar fields on static spherically symmetric spacetimes. 

We now consider the electromagnetic analog of the above calculations. Since the calcu- 
lation is highly analogous, we simply sketch the procedure and present the results. For a 
static worldline z(t) on our spacetime (jlj), the only contribution of the Green's function to 
the self-force (130]) is from the time-time component, which is scalar under spatial rotations. 
We may therefore expand in ordinary scalar harmonics, as in the case of a scalar charge 
treated in section IIHI The treatment then proceeds in precise analogy, and we find that for 
an electric charge held fixed outside an arbitrary (nonpolarizable) body of the form (jlj), the 
self-force difference from pure Schwarzschild is given by 

J2 



with 



1=0 



^ W(I e ,P e 



W(I t ,Q e ) 



(38) 



r=r 

and where (as always) derivatives in the Wronskian are taken with respect to r (rather than 
f or f). Here Pg and Qe are the electrostatic mode functions for Schwarzschild spacetime 
(i.e., the radial parts of the pure multipole solutions to the electrostatic Maxwell equation 
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on Schwarzschild spacetime), given by derivatives of Legendre functions 33], 



I I(Th) ( r ~ 1 ) P £ ( r ) otherwise 
Q i (f) = (f-l)Q' e (f), 




(39) 



(40) 



and li is the origin-regular electrostatic mode function for the interior spacetime, given by 
the (unique up to normalization) origin-regular solution of 




Equations fl37j) . fl38|) . and (TJT]) are the electromagnetic analogs of equations (135]) . fl34|) . and 
(IE]) (specialized to to — 0), respectively. Note the change in sign of the \&' term in equation 
( 14T1) from its scalar analog (|H]). Note also the lack of Ricci terms in equation f j4TT) (and 
corresponding lack of a boundary curvature term in equation fl38l ). despite the presence of 
Ricci terms in the equation (I3T1) . 

In the far- field limit r q — » oo, again the £ = mode dominates the self-force difference on 
account of the r~^ +1 ** behavior of Qe- However, as in the minimally-coupled scalar case, both 
Pq and Jo are just constants, so that the Wronskian W(Io, Po) vanishes and the coefficient 
of the leading order term is in fact zero. Thus the self- force difference vanishes at leading 
order, and we have shown that the self-force on a static electric charge is independent of the 
central body type at leading order in the far field, confirming the suggestion of table HI 

If we restrict momentarily to smooth central bodies (i.e. those with no boundary layer 
at r = To), then in each case investigated so far we have found that the leading-order far- 
field static self-force difference is given by the Wronskian of the u = i = origin-regular 
interior solution with the u = £ = horizon-regular Schwarzschild solution. Since the 
calculations done here in the scalar and electromagnetic cases generalize straightforwardly 
to other separable wave equations, it is clear that this conclusion must apply quite generally. 
When asking if the far- field static self-force will be independent of the central body for a given 
wave equation, therefore, one may ask if the regular static spherically symmetric solution for 
any central body always agrees with that of the Schwarzschild metric. For example, in the 
case of a massive scalar field, it is easy to see that the origin-regular spherically symmetric 
solution will be different in different spherically symmetric spacetimes, so that weak locality 
will not hold in this limit. This is clearly the general case; the minimally-coupled scalar 
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and electromagnetic wave equations have the quite special property that the regular static 
spherically symmetric solution is a constant for all metrics. Therefore, it seems natural to 
view the weak locality of the static far-field force as a special property of the minimally- 
coupled scalar and electromagnetic cases[j^ 

IV. THIN-SHELL SPACETIME: STATIC AND CIRCULAR ORBITS 

We now adopt the specific choice of a central body composed only of a thin shell of 
matter, which enables the calculation of concrete results. That is, we choose the interior 
spacetime to be flat, so that the full metric is 

{-f(r )dt 2 + df 2 + f 2 dVt 2 < f < r 
(42) 
-fdt 2 + j- l dr 2 + rW r > r 

In the language of equation we have chosen e 2 * = f(r ) and e 2A = 1. For these choices, 
equation (j5J) becomes 

f =r + (r-r )/(r )- 1/2 . (43) 

It is easy to see that while the metric components in the coordinates t, r, 8, <ft are continuous 
functions across the boundary r = tq, their r-derivatives are not. Their second r-derivatives 
therefore contain a delta function, which is interpreted as the spacetime curvature due to 
a thin shell of matter present at r = ro jll, [22}. Our interest is in the Ricci scalar of the 
spacetime, which we calculate to be 

R = X 6(r - r ) (44) 

with 



2 

X = — 

r 



^2v7(^-2/(r )-^V (45) 

The quantity T = — ( y l/8n)f~ 1 / 2 x has the interpretation of the trace of the surface stress- 
energy tensor of the shell, with the factor of f^ 1 ^ 2 correcting for the fact that x is the 
coefficient of a delta function in r, rather than in proper distance. For ro — > 00 we have 
— T — > (l/47r)M/ro, giving the surface density of a Newtonain shell. For r — > 2M, we 

14 We note that the gravitational self-force does not make sense in the static case, since no non-gravitating 
strut is available to hold a mass in place. If one simply computes the field of a static charge, this manifests 



itself as a conical singularity in the solution 



34|. 
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have T — > oo, showing that there cannot be a shell as small as its Schwarzschild radius. 
Finally, at ro = (9/4)M we have T = 0. It is curious that this occurs at ro = (9/4) M, 



which is the minimum size established by Buchdahl [36] for smooth interiors. (For example, 
a constant density star develops a pressure singularity for r /M = 9/4.) For our purposes 
the relevance of r = (9/4) M is as 1) the size of shell where the force is independent of the 
coupling to curvature and 2) the maximum size for which instabilities can occur for £ > 
(see appendix). We note that our result for the Ricci scalar agrees with the (more general) 



results of 



35]. 



A major advantage of the thin-shell spacetime is that — since the interior metric is flat- 
its mode functions are known in closed form, even in the non-static case. Equation (JSJ) for 
the mode functions becomes simply 

G'Ur) + ?GUf) - (fcil - /(r o rv) GUr) = 0, (46) 

which is of course the ordinary scalar wave equation for flat spacetime. The solution regular 
at the origin is simply (e.g. |23| ) 



r l u = 

lUr) = { , (47) 

J e (urfiro)- 1 ^ 2 ) u^O 

where Je(x) gives the spherical Bessel function regular at x — 0. 
A. Static Case 

With exact solutions for the interior mode functions at hand, we can evaluate fl35l) ex- 
plicitly to determine the self-force on a static charge. We reproduce the results of [30J for 
£ = 0, and find that £ > increases the repulsive self-force. (We exclude £ < 0, which can 
give rise to instabilities (see appendix) even for large shells. We also exclude r < 2.25M 
for £ > to avoid instabilities.) We note that the self- force diverges as the particle ap- 

n 

proaches the shell, as found for minimal coupling in [30J. Our main interest in the static 
case, however, is in the far-field self-force (136|) . Using the explicit forms Iq = 1,Pq = 1, and 
Qo(x) = l/21og[(x + l)/(x — 1)], the leading order piece of the self-force becomes 

F far- fie id = ^r 2 oX (r ) (4g) 

r l 2M + £r 2 x(r ) log / 



t r -2M 
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Equation ( HHj) gives the self-force on a static charge held a large distance away from a 
spherical shell of radius r$. Two limits are of interest. First, when the shell radius tq 
becomes very large, the results should go over to those of Pfenning and Poisson [14( for 
Newtonian bodies. Indeed, it is easy to see that we have 

F far-ficid ^ 2 ^y_ ag r0 ^ OO; (49) 

in agreement with their results. Second, when the shell radius tq is very nearly 2M, the 
result should approach Wiseman's result [16| of zero self-force in the black hole case. Indeed, 
it is easy to see that we have 

F far-fieid^ ag rQ ^2M. (50) 

Therefore, our results do reproduce the known results in the limiting cases of weakly curved 
(Newtonian) and highly curved (black hole) central bodies. To see the regime intermediate 
between these extremes, we can plot equation fHHj) directly (figure [2]). This illustrates the 
change in far-field self-force as the spacetime transitions from one that is weakly curved 
everywhere to a spacetime with a highly curved central object. For £ > 0, the self-force 
is repulsive for most of parameter space, becoming attractive only in the dubious regime 
r < 2.25M where instabilities and arbitrarily large self-forces can occur (see appendix). The 
sensitivity of the far-field self-force to the central object type dramatically demonstrates the 
non-locality of the tail integral. 



B. Circular Orbits 



The four-velocity of a circular (equatorial) orbit of radius r q and frequency Q is given by 
u a = it*(l,0, 0, fl) with u l = (f(r q ) — Vt 2 r 2 )~ 1 / 2 . For geodesic orbits (which we consider), 
one has the radius-frequency relationship Q = M/r^ allowing one to simplify this to 
it* = 1/ a/1 — 3M/r q . The worldline coordinates are z M (r) = (wV, r q , tt/2, fhtV). Evaluating 
equation fl3]) given the decomposition (J7J) then gives 



AF M = — 



E ^(ro)iV to i2S(r)^(r 3 )ZrPr(cos0)Pf(O) e ^- nt ) 

=0 m=-£ 



(51) 



x n =z n 
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FIG. 2. A plot of the coefficient of the far-field force, equation (|48|). for £ = 1 as the shell radius 
ro is varied. For very large r$ (not shown), the coefficient asymptotes to 2, the Newtonian result. 
For ro close to 2M, the coefficient approaches the black hole result of zero (albeit with infinite 
slope). The force changes sign at ro = 2.25M, and becomes arbitrarily large near an instability 
(see appendix) present for ro ~ 2.09M. 



where we have 

vm (2l + l)(£-m)\ 



(52) 



1 Air (£ + m)\ 

and Cj = mVt. The integral over frequencies collapsed to a discrete sum on account of the 
harmonic dependence on the t and <fi cordinates. Component-by-component, equation ( I5T1) 
gives 

A F t = -VlAF^ (53) 

2 00 i 

AF r = 77 E E WWW ^ W (0) 2 ^(ro) (54) 



£=0 m=-£ 

AF e = (55) 

2 oo £ 



AF ^ = ^tJ2Yl i™N^BZ{r q ) 2 Z?Pr{Q) 2 KUr )- (56) 

£=0 m=-£ 

Note that the summands in equations fl54|) and fl56l) complex conjugate under m — >■ — m, 
implying that only their real parts contribute to the overall (real) sum. 

Equations f )53H56|) give the self-force difference from pure Schwarzschild for a particle in 
circular orbit about an arbitrary central body of the form fl3j). For the thin-shell spacetime, 
the interior mode functions I(& are given by ( 14T|) . so that the only quantities not known 
in closed form (or readily available numerically) are then the Schwarzschild mode functions 
and Rf m . We therefore compute these numerically, by solving equation fl9]) subject 
to the boundary conditions f lT6|) and ( ITT)) . One uses these equations to provide initial 

22 



valued_| far from the particle, and then numerically integrates equation fl9]) to the required 
r q . We used the software package Mathematica 7 [37J for all of our numerical computations. 
We have verified the accuracy of our mode functions by using them to compute the fluxes 
reported in 39] (which have now been independently confirmed in 40j) in the special case of 
Schwarzschild. We have also verified that the flux agrees with the local dissipative self-force 
for all radii, in the manner described in (e.g.) 39]. 

With all the ingredients assembled it is now a simple matter to evaluate equations f )54|) 
and ( 156]) for various values of £,,r ,r q . As in the static case, we exclude £ < as well as 
£ > 0,r < 2.25M to avoid instabilities. In general the self-force difference (radial and 
angular) is positive, increasing the radial force found in 4l|, and decreasing the magnitude 
of the negative angular force (e.g. [39|). As £ is increased (at fixed r ,r q ), the (radial 
and angular) self-force difference increases. Likewise, as r$ is increased (at fixed r q ,£), the 
(radial and angular) self-force difference increases. As in the static case, the radial self- 
force diverges as tq approaches r q (while the angular self-force remains finite). As r q is 
increased (at fixed £,ro), the magnitude of the self- force difference decreases. Since we find 
no qualitative differences for £ > (and ro > 2.25M), for simplicity we will restrict to 
minimal coupling (£ = 0) for the remainder of this section. 

Figure [3] lists the results for various values of r and r q , showing some of the trends 
discussed above. In order to better understand the results, however, it is instructive to 
instead consider the fractional self-force difference; that is, the magnitude of the self-force 
difference divided by the self-force in one of the spacetimes. The self-force on a scalar charge 
in circular orbit about a Schwarzschild black hole was computed in 411] • The authors list 
only the radial self-force, but the angular self-force may be easily computed from the mode 
functions Rf w and R\ w as described (e.g.) i n |39 ]. (The radial self-force, on the other hand, 
requires regularization.) Using data from 41( for the radial black hole self- force and our 
own computations for the angular black hole self-force, we compute \AF r / F^ lack hole | (and 
likewise for 0) for a variety of radii (figure HI). The main advantage of the fractional self- force 
difference is that the radial and angular components can be sensibly compared. We see that 
the radial fractional self-force difference is always order unity, while the angular fractional 
self-force difference is small and falls off rapidly. 

In fact, we use a higher-order series representation to provide the initial values, in order to allow the 
integration to begin closer to the particle. Higher-order terms in the end behavior (1161) and (1171) are 



computed from a recursion relation that follows from equation ©, as described (e.g.) in [38[ and 
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FIG. 3. A table of results for the self- force difference, equations (j54|) and (|56p . in units where 
M = 1 and with £ = and q = 1. The first /last three rows show the radial/angular self- force 
difference for three different shell radii r$. The columns vary the particle radius r q . 
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FIG. 4. The fractional self- force difference \AF/F 



black hole 



, for radial (conservative) and angular 



4l|.) Here we have used 



(dissipative) components. (Black hole self-force results are taken from 
units where M = 1 and taken £ = 0, ro = 5. The change in central object is seen to have a much 
larger effect on the radial force than on the angular force. The increase in the radial component 
near the shell is associated with the divergence of the radial shell self- force when r q = tq. 



This result can be understood by recalling that the radial force is conservative, while the 
angular force is dissipative. More specifically, with G re t/adv the retarded/advanced Green's 
function, one can define Gdiss = (l/2)(G rc t ~ C a d v ) and G cons = (l/2)(G rct + G a dv), so that 
we have G ret = Gdiss + G cons . A self- force computed with G CO n S will be time-symmetric and 
therefore cannot carry energy away from the system (hence the name conservative); thus 
the remainder is entirely responsible for energy loss (hence the name dissipative). Noting 
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that the interchange ret-B-adv corresponds for circular orbits to t ■<-> —t, <fi -H- —0, it is 
easily established that the radial force is purely conservative, while the angular force is 
purely dissipative. This puts the results of figure H] in context: It makes sense that the 
dissipative piece would be dominated by local emission of radiation, with only small effects 
due to the excitation of the central bodyj 16 ! On the other hand, there is no reason why the 
conservative piece should not depend strongly on the central body, especially in light of the 
fact that conservative self-forces often owe their entire existence to boundary conditions. For 
example, the force on an electric charge held static outside a conductor would be counted 
as a (conservative) self-force under our definitions. 

The result of figure H] — combined with the above interpretation — suggests that in general 
the dissipative part of the self- force should be more weakly local (i.e., independent of the 
central body type) than the conservative part, even in the strong field regime. It would be 
very interesting if the dissipative self-force could also be shown to be more truly local (in the 
sense of tail falloff) than the conservative part. While this cannot be true in the weak field 

nn 

given the computations of |13|, H4j , it seems possible that it could become true for near-field 
orbits. In general, it seems plausible that the tail could become more local in the near-field, 
where the light reflection time is decreased. 

An interesting question to ask is to what extent the circular orbit conservative self-force 
agrees with the static case self-force. In pure Schwarzschild spacetime there is no agreement: 
the static case self-force vanishes, while the circular orbit force does not. Likewise, the static 
and circular orbit self-forces do not agree for the thin-shell spacetime, differing by factors 
of order unity. However, the self-force difference does in fact closely agree, as seen in 



figure figure 
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r q = 20M 



5l The good agreement even at reasonably small radii (at the percent level by 
indicates that the circular orbit radial self-force in the shell spacetime may be 
estimated by taking the circular orbit results in pure Schwarzschild spacetime, and adding 
to them the self-force (difference) computed for the shell in the static case. Combined with 
the result that the dissipative self-force difference is very small, this suggests that (except 
in the very strong field) only static case results are necessary to convert black hole results 
to results for other bodies. 
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Indeed, one can check that the Newtonian result of 



14] which explicitly takes a radiation- reaction 



form — gives the correct dissipative self-force (for either spacetime) at the percent-level for r q = 30M, at 

the five-percent-level for r q = 15M, and still within fifty percent at r q = 6M. 
17 If £ or ro is increased, the radius of agreement moves outward somewhat. 
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FIG. 5. A plot of the ratio of the circular orbit radial self-force to the static case radial self-force, 
along with the same quantity for the self-force difference (which agrees with the self-force in the 
static case, Ai^ tatlc = F^ tatlc ). We see that there is disagreement in the self- force even for large r q , 
while there is very good agreement for the self-force difference. The results of [41] have been used 
in the calculation of these quantities. We have used units where M = 1, and taken tq = 5, £ = 0. 



V. SUMMARY 



For a particle orbiting a static spherically symmetric body, we have studied the depen- 
dence of the self-force on the choice of interior metric. We began with a general treatment of 
the problem in the scalar case (additionally outlining the electromagnetic and gravitational 
cases), and then considered specific central bodies and/or orbits. We first considered static 
charges in the scalar and electromagnetic cases, confirming the suggestion of earlier work 
[l^ that the far-field self-force is independent of the central body type in the minimally 
coupled scalar and electromagnetic cases, but dependent on the central body type in the 
nonminimally coupled scalar case. We then adopted the specific choice of a thin-shell central 
body and computed the self-force difference from a black hole central body in the case of 
static and circular orbits of scalar charges. We found that the change in central body has 
a much larger effect on the radial (conservative) self-force than on the angular (dissipative) 
self-force difference. We also found that the radial self-force difference is well approximated 
by the static case (radial) self-force difference, raising the possibility of using static case 
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results to correct self-forces (for arbitrary orbits) for a change in central body. 
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Appendix A: Instabilities 

If the denominator in the expression fl23|) for the self-force difference vanishes, then that 
expression becomes ill-defined, indicating (as remarked below (I2T]) ) that the mode decom- 
position ([7]) is not defined. We now show that the origin of this problem — if it occurs — is 
a linear instability of the offending spcetime. Since the denominator is proportional to 
of equation ( 1281) . it is clear by equation ( 1271) that it vanishes only for the special case where 
the interior solution Ii w matches directly to the infinity solution Kf° (with zero coefficient of 
the horizon solution Rf u ). This would give a homogeneous solution to scalar wave equation 
on the spacetime (TjJ that is bounded and purely outgoing at infinity. If u ^ 0, it is easy 
to check that any such solution would violate conservation of total energy, and therefore 
cannot exist. However, we can have such behavior when u = 0, which corresponds to the 
existence of a static homogeneous solution that is bounded. In this case multiplication by t 
gives a second spatially-bounded homogeneous solution, so that the spacetime has a linear 
instability. While it may be possible to make sense of a self-force in such a spacetime, we 
make no attempt in this paper. 

However, it is still useful to categorize — as best as possible — when this type of instability 
can occur. This is a straightforward matter for the I = mode of the thin-shell spacetime, 
for which explicit expressions were worked out in section IIVt and we carry this out below. 
The main conclusions are that the product £x must be negative (corresponding to a negative 
mass-squared term in the Klein-Gordon equation) for an instability to occur; and for £ > 0, 
this can happen only for shells in the extreme region tq < 2.25M. We are unable to provide 
analogous results for £ > 0, although it would seem reasonable to expect the same properties 
to hold. 
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We will use the properties of x( r o) (equation f l4"5]) ) that 1) x is bounded on (2M, oo); 2) 
X has one zero at r$ = 2.25M, being negative on (2M, 2.25M) and positive on (2.25M, oo); 
and 3) x bas one local maximum at 3M, being increasing on (2M, 3M) and decreasing on 
(3M, oo). Now, the denominator of equation ( l4"gj) vanishes — and hence there is an I = 
instability — if and only if 



Thus, C,x must be negative for the instability to occur. 

To analyze in more detail, consider first the case of £ > 0. In this case h(r ) is negative, 
monotonic, and unbounded on (2M, 2.25M). This implies that it intersects —1 exactly once 
in that domain. For r$ > 2.25M, on the other hand, h(ro) is positive and does not intersect 
— 1. Therefore, for £ > there is exactly one radius for which an instability occurs, and this 
radius is bounded by 2M < r < 2.25M. 

Now consider £ < 0. This flips the sign of h(ro) so that it may intersect —1 only for 
r > 2.25M. Since h(ro) is bounded and contains a single local maximum (at r = 3M), this 
means that there may be zero, one, or two intersections with —1. The value of £ for which 
there is one intersection is determined by h(3M) = —1, giving £ »s —.98. For —.98 < £ < 
there is no intersection, and for £ < —.98 there are two intersections. 



[1] E. Poisson, Liv. Rev. Rel. 7 6 (2004) 

[2] B.S. DeWitt and R.W. Brehme, Ann. Phys. (N.Y.), 9, 220-259, (1960) 
[3] T.C. Quinn, Phys. Rev. D 62 064029 (2000) 

[4] S.E. Gralla and R.M. Wald, Class. Quantum Grav. 25 205009 (2008). 
[5] S.E. Gralla, A.I. Harte and R.M. Wald Phys. Rev. D 80, 024031 (2009) 
[6] A.I. Harte, Class. Quant. Grav. 25 235020 (2008) 
[7] S.E. Gralla, Phys. Rev. D 81 084060 (2010) 
[8] E. Poisson, |arXiv:gr-qc/9912045| (1999) 

[9] P. R. Anderson and B. L. Hu, Phys. Rev. D 69, 064039 (2004) 
[10] W.G. Anderson, E.E. Flanagan, A.C. Ottewill, Phys. Rev. D 71, 024036 (2005) 
[11] M. Cassals, S. Dolan, A.C. Ottewill and B. Wardell Phys. Rev. D 79, 124043 (2009) 
[12] W. G. Anderson and A. G. Wiseman, Class. Quant. Grav. 22, S783 (2005) 




(Al) 



28 



[13] B.S. DeWitt and CM. DeWitt, Physics (Long Island City, N.Y.) 1, 3 (1964) 

[14] M.J. Pfenning and E. Poisson Phys. Rev. D 65, 084001 (2002) 

[15] A. G. Smith and C. M. Will, Phys. Rev. D 22, 1276 (1980) 

[16] A.G. Wiseman, Phys. Rev. D 61 084014 (2000) 



[17] [http://lisa.nasa.gov 



[18] L. Barack, Class. Quant. Grav. 26 213001 (2009) 

[19] R.M. Wald 1984 General Relativity (Chicago, IL: University of Chicago Press) 

[20] B.F. Schutz 1985 A first course in general relativity (Cambridge, UK: University Press) 

[21] E. Poisson 2004 A Relativist's Toolkit (Cambridge, UK: University Press) 

[22] W. Israel, Nuovo Cim. B 44 1 (1966) 

[23] J.D. Jackson 1999 Classical Electrodynamics (Hoboken, NJ: John Wiley & Sons) 
[24] K.D. Kokkotas and B. Schmidt, Living Rev. Relativity 2 2 (1999) 

[25] W.E. Boyce and R.C. DiPrima 2004 Elementary Differential Equations and Boundary Value 

Problems (New York, NY: John Wiley & Sons) 
[26] T.C. Quinn and R.M. Wald, Phys. Rev. D, 56, 3381-3394, (1997) 
[27] Y. Mino, M. Sasaki, and T. Tanaka, Phys. Rev. D 55 3457 (1997) 
[28] A. Pound, Phys. Rev. D 81 024023 (2010) 
[29] K. Shankar and B. Whiting, Phys. Rev. D 76 124027 (2007) 
[30] LM. Burko, Y.T. Liu and Y. Soen, Phys. Rev. D 63, 024015 (2000) 
[31] G.B. Arfken and H.J. Weber 2001 (San Diego, CA: Academic Press) 
[32] P. Anderson and A. Fabbri, Phys. Rev. D 75 044015 (2007) 
[33] J.M. Cohen and R.M. Wald. J. Math. Phys. 12 1845 (1971) 
[34] T.S. Keidl, J.L. Friedman and A.G. Wiseman Phys. Rev. D 75 124009 (2007) 
[35] A. Eid, Electronic Journal of Theoretical Physics, 5, No. 19 115-122 (2008) 
[36] H.A Buchdahl Phys. Rev. 116 1027 (1959) 

[37] Wolfram Research, Inc., Mathematica Version 7.0, Champaign, IL (2008) 

[38] S. Detweiler, E. Messaritaki and B.F. Whiting Phys. Rev. D 67, 104016 (2003) 

[39] S.E. Gralla, J.L. Friedman and A.G. Wiseman, [arXiv:gr-qc/0502123 



[40] N. Warburton and L. Barack Phys. Rev. D 81, 084039 (2010) 

[41] L.M. Diaz-Rivera, E. Messaritaki, B.F. Whiting and S. Detweiler Phys. Rev. D 70 124018 
(2004) 

29 



